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Since the discovery of adiabatic quantum computing, a need has arisen for rigorously proven 
bounds for the error in the adiabatic approximation. We present in this paper, a rigorous and 
elementary derivation of upper and lower bounds on the error incurred from using the adiabatic 
approximation for quantum systems. Our bounds are often asymptotically tight in the limit of slow 
evolution for fixed Hamiltonians, and are used to provide sufficient conditions for the application 
of the adiabatic approximation. We show that our sufficiency criteria exclude the Marzlin-Sanders 
counterexample from the class of Hamiltonians that obey the adiabatic approximation. Finally, we 
demonstrate the existence of classes of Hamiltonians that resemble the Marzlin-Sanders counterex- 
ample Hamiltonian, but also obey the adiabatic approximation. 



I. INTRODUCTION 



The adiabatic approximation is a technique of central importance to the study of time-dependent quantum systems. 
Its importance is derived from the fact that adiabatic evolution is one of the few known limiting cases where the 
. dynamics that is generated by a time-dependent Hamiltonian is well understood. The approximation states that 
if the Hamiltonian that generates a quantum system's evolution changes sufficiently slowly, then a quantum state 
I""!, that is prepared in an instantaneous eigenstate remains in an instantaneous eigenstate throughout the evolution. 

■ The adiabatic approximation has been used to control the evolution of quantum systems [lillli analyze particle 
collisions Q , design circuit based quantum computers and has been used to create a new paradigm for quantum 

tl^ computing; namely adiabatic quantum computing 0, |J 
^ In many of these applications, bounding the time needed to apply the adiabatic approximation was seen as a 

^ 1 secondary issue. Adiabatic quantum computing proved to be a notable exception. In adiabatic quantum computing, 
O^' quantum algorithms are constructed by choosing a time-dependent Hamiltonian whose ground state at t = is an 
easily prepared quantum state and at < = T, the ground state encodes the solution to a computational problem. If 
the adiabatic approximation holds, then the time-evolution transforms the initial ground state into the solution to the 
I problem, thereby enacting the quantum algorithm. A different notion of algorithmic complexity must be used to assess 

■ these algorithms because they do not use a discrete set of gates. The scaling of the energy and the time needed for 
' this evolution with the problem size, is a natural way to assess the complexity of such algorithms jzl, i9|]. Consequently, 

rigorous estimates for the error in the adiabatic approximation are needed to demonstrate such a scaling [lo| . 

There are many possible ways to adiabatically evolve the initial state into the final state [lll - fTsj . We are, at least 
, in principle, given the freedom to choose these Hamiltonians to depend on the properties of the evolution such as 
' the Hilbert space dimension, the error tolerance or even the total evolution time. This line of inquiry lead to the 
_ Marzlin-Sanders counterexample [l4ll , which overturned conventional wisdom about the validity of the approximation. 
I • In 2004 Marzlin and Sanders [3] put forward a Hamiltonian that satisfied the traditionally held "slowness" crite- 
^ ' rion for the validity of the approximation, and yet it generated an evolution that did not agree with the adiabatic 
approximation, even in the limit of arbitrarily slow evolution. This created doubt in the scientific community about 
the conditions that are necessary for the adiabatic approximation and, along with the discovery of adiabatic quantum 
, computing, has led to further interest into the conditions are needed to apply the adiabatic approximation |15l4l7l |. 
We achieve three important goals in this paper. First, we provide a rigorous analysis of adiabatic evolution that is 
intuitive, elementary and avoids ad hoc physical assumptions. Second, we provide upper bounds and lower bounds 
for the error in the adiabatic approximation that are often asymptotically tight, in the limit of slow evolution for a 
fixed Hamiltonian. Third and last, we show that the adiabatic approximation can be successfully applied to a broad 
class of quantum systems that are related to the Marzlin-Sanders counterexample. 



II. BACKGROUND 



In this section we present the background material that is needed to understand our main results, which are given 
in the following section. In particular we discuss the adiabatic approximation in greater detail, the Marzlin-Sanders 
counterexample 14 1 and the Jansen, Ruskai and Seller (JRS) bound on the error in the adiabatic approximation [l0|. 

The time-evolution of closed quantum systems over a time interval t e [0,r] is described by the time-evolution 
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operator, which is expressed in terms of the dimensionless time s = t/T as 

u{sMm)--=ms)), (1) 

where \^{s)) is a quantum state evaluated at the dimensionless time s. 

We write the adiabatic approximation for the action of the time-evolution operator on an instantaneous eigenstate 
of H{0), |G(0)) (typically taken to be the ground state) as 

C/(s,0)|G(0)) « e-'^°^««)^«^|G(s)), (2) 

where Eg{s) is the eigenvalue of the instantaneous eigenvector |G(s)), and the geometric phase [l^ is absorbed into 
the definition of the eigenstate. The error in the adiabatic approximation is the magnitude of the projection of the 
evolved state onto the space that is orthogonal to |G(s)). This error is generally considered to be small if the variation 
of the Hamiltonian with time is small compared to the square of the minimum eigenvalue gap between the initial 
state |G(s)), and all other instantaneous eigenvectors [l^,[23]- We write this requirement as 

max,|lij(.)|| ^^^^^^ 



mills, i,[j„,Gis)] 

where ju,Gis) is the eigenvalue gap between the instantaneous eigenstates \iy{s)) and |G(s)). Adiabatic theorems have 
also been proven for systems that have eigenvalue crossings, but approach those crossings in a smooth fashion f2l| . 

The Marzlin-Sanders counterexample [14| provided a Hamiltonian for which Eq. [3] holds but Eq. [2]does not. Their 
counterexample Hamiltonian is chosen such that ||-ff(s)|| diverges in the limit as the quantum system's evolution 
becomes arbitrarily slow. The discrepancy with Eq. [5] was attributed to a resonant effect, wherein the normally 
negligible contribution of the second derivative of H became dominant due to increasingly rapid variation as T 
increases. Although the fact that the error in the adiabatic approximation depends on the differentiability class of the 
Hamiltonian was recognized before the Marzlin-Sanders counterexample was discovered [22rf24| . the counterexample 
showed that terms that depend on higher-order derivatives of H{s) can sometimes dominate the expression for the 
error in the adiabatic approximation. 

This problem can be addressed by using upper bounds for the error in the adiabatic approximation, rather than 
asymptotic estimates, because upper bounds do not neglect any contribution to the error. Examples of such upper 
bounds can be found in (lol.[2 ^ -|25| . The best upper bound presently known for the error in the adiabatic approximation 
is given by the JRS bound [10] ■ The JRS bound states that for an to- fold degenerate initial state |G(0)) (meaning that 
there exist to — 1 instantaneous eigenvectors that share the same eigenvalue as |G(s)) for all s), the approximation 
error ||(1 - Fg(1))C^(1, 0)|G(0)) || is bounded above by 



m||ff(0)|| m||ij(l)|| 1 ( TO,||H(0|| 7my/^\\H{^)\\ 
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de (4) 



The quantity 7(^) in Eq. |4] represents the minimum eigenvalue gap between the state |G(0) and any other instan- 
taneous eigenstate at time s = ^, Pg{^) represents the projector |G(1))(G(1)| and 1 is the identity operator. Given 
that the Hamiltonian -ff(s), is not explicitly a function of T, this upper bound can be inverted to find a value of T 
such that the error in the adiabatic approximation is less than any fixed error tolerance. 

As Eq.|4]depends on ||if(s)||, it seems to be ideal for addressing the issues caused by the divergent second derivative 
of the Marzlin-Sanders counterexample Hamiltonian. However, JRS noted that the third term in Eq.|3]can be replaced 
by one that is 0(1/T^) [l^illll, implying that their upper bound is not asymptotically tight. Because their bound is 
not asymptotically tight, it cannot be used to find a sufficient condition for the validity of Eq. [S] 

We expand on these results by providing upper and lower bounds for the error in the adiabatic approximation 
that are often asymptotically tight (meaning that our error bounds converge to the observed error in the limit of 
large T), given that the leading order term in our bounds does not vanish for any T. We then use these bounds 
to address Hamiltonians that are related to the Marzlin-Sanders counterexample. We find using these bounds that 
many Hamiltonians whose second derivatives diverge as T increases, obey the adiabatic approximation. These results 
are summarized in the following section. 



III. RESULTS 



This section contains our main results, which address the question of when the adiabatic approximation can be 
used and provide asymptotic estimates of the error in the approximation. These results are consequences of our upper 
bounds for the error in the adiabatic approximation, which we present in Section [VHI 
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As with most adiabatic theorems, we require that the variation of the Hamiltonian is slow compared to an appro- 
priate timescale. Our first theorem states that the adiabatic approximation is vahd if the evolution time T, is large 
compared to a timescale Aq, which we define below. Our timescale Ag reduces to that in Eq. |3]if T is large compared 
to a second timescale Ai, which depends on the second and third derivatives of the Hamiltonian. These timescales 
are given in the following definition. 

Definition 1. Let H{s) be differentiable three times on the interval [0,1], ^v,Gis) be the eigenvalue gap between 
the instantaneous eigenstates \G{s)) and \v{s)) and let jmin be the minimum eigenvalue gap between any two non- 
degenerate instantaneous eigenvectors of H{s) on this interval. We then define the timescales Aq and Ai to be 



Ai := 



\H\\ 



Ai 



7n 



^3 . ^6 ^ 

' mm ' mill 



(5) 
(6) 



where \\H\\, \\H\\ and \\H\\ represent the maximum values of the norms of the first three derivatives of H[s). 

The following theorem uses the criterion that the timescale Aq is asymptotically smaller than T as a sufficient 
condition for the applicability of the adiabatic approximation, in the limit of large T. We refer to the resulting 
theorem as a zeroth-order adiabatic theorem because it gives an expression for the approximation error that is correct 
to zeroth-order in powers of . 

Theorem 1. (Zeroth-order adiabatic theorem) Let H : [0,1] i— > C^^^ be a Hamiltonian that is differentiable three 
times and has a minimum eigenvalue gap between any two non- degenerate eigenvectors of jmin > 0. // \G{s)) is a 
non-degenerate instantaneous eigenstate of H{s) and Ag G o{T), then the error in the adiabatic approximation obeys 



{1-Pg{1))U{1,0)\G{0))^0{Ao/T), 



(7) 



where 1 is the identity operator. 



The above theorem states that the error in the adiabatic approximation is zero, plus an additional term that is 
small in the limit of large T. The requirement that Ag G o{T), reduces to the standard criterion given in Eq. [3]if the 
derivatives of L[{s) are independent of T. This implies that Theorem [1] is as a generalization of the standard criterion 
for the validity of the approximation. The following theorem improves upon the result of Theorem [1] by providing an 
expression for the error in the adiabatic approximation that is correct to leading order in powers of T~^. We call this 
result a first-order adiabatic theorem. 

Theorem 2. (First-order adiabatic theorem) Let H : [0, 1] i— )■ C^^^ be a Hamiltonian that is differentiable three times 
and has a minimum eigenvalue gap between any two non- degenerate eigenvectors of"f^in > 0. If we take \G{s)) to be a 
non-degenerate instantaneous eigenstate of H(s) and Ai g o(7niin2^^); then the error in the adiabatic approximation, 
(i/-PG(l))C/(l,0)|G(0)),, IS 



-«7G,i/(si)r 



o 



Ai 



y2 



(8) 



si=0 



where = E^{£_)d£^T, ^g,u{s) Eg{s) — Eu{s), -Pg(I) = IG'(l)) (G(l)| and the phase of each instantaneous 
eigenstate is chosen such that {v{s)\i'{s)) = for every s G [0, 1] and every v G {0, . . . , — 1}. 

It may seem strange that the leading order expression in Eq. [5] only depends on the properties of the Hamiltonian 
at the beginning and end of the adiabatic evolution, in contrast to Eq. [31 This dependence occurs because the time- 
evolution operator is expressed as a sum of rapidly oscillating path-integrals in the adiabatic limit [1^. We expect 
that the contribution of each complete oscillation will have a negligible contribution to the integral if the oscillation 
frequency is sufficiently rapid compared to the timescale for the variation of the Hamiltonian. Therefore we expect 
that the contribution of the incomplete oscillations will dominate the behavior of the integral. Only the first and last 
such oscillations are not necessarily complete. The contribution of the first the last and last oscillation is dominated 
by the properties of the Hamiltonian at the beginning and the end of the evolution respectively; therefore, we expect 
that the leading order expression for the error in the adiabatic approximation should be dictated by the pro perties of 
the Hamiltonian at those times. This dependence is also noted by a number of other studies (Tol. [isl [26l l27j]. 
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These theorems imply that the adiabatic approximation also applies for some Hamiltonians that behave like the 
Marzlin-Sanders comiterexample in that ||-ff|| and \\h\\ are increasing functions of T. Our first corollary, which is 
given below, provides a criterion for the size of the derivatives of H that is sufficient to guarantee that the standard 
estimate of the error in the adiabatic approximation applies. 

Corollary 1. If H : [0,1] i-^ ([^nxn dijferentiable three times, \G{s)) is a non- degenerate eigenstate of H{s), 
\\H\\ G o{^/T), G o{T) and both \\H\\ and (7min)~"'^ o,i"e bounded above by a constant function of T for all T > 
then 

11(1 - Fg(1))C/(1, 0)|G(0))|| e0( max 

\ 5=0,1 

In other cases the adiabatic approximation may hold, but the error in the adiabatic approximation may not be 
0{1/T). A sufficient condition for the validity of the adiabatic approximation, for a sufficiently large T, is given 
below. 

Corollary 2. If H : [0, 1] i— > C^^^ is twice dijferentiable, \G{s)) is a non-degenerate eigenstate of II{s), \\II\\ G o{T) 
and both \\H{s)\\ and (7min)^^ are bounded above by a constant function of T for all T > then 



\\His)\\ 
niini.[7j.,G(s)]^T 



(9) 



lim |1(1-Pg(1))C/(1,0)|G(0))|| =0. (10) 

T— >-oo 

A principal limitation of these results is that we take \G{s)) to be non-degenerate for all s. We make this assumption 
to simplify our analysis of the time-evolution operator, and leave the general case for future work. 

These results are laid out as follows: in Sec. II VI we discuss the concept of paths in greater detail and their relation 
to time-evolution. In Sections |V] and IVII we use the results of Sec. IIVI to provide a path integral representation of 
the time-evolution operator. In Sec. IVIII we use integration by parts to approximate the path integrals and prove 
Theorems [T] and [21 We then compare our error bounds to the error found for the Search Hamiltonian in Sec. I Villi 
and discuss the Marzhn-Sanders counterexample [U] and prove Corollaries [5] and [1] in Sec. IIXI 



IV. PATHS 

The concept of paths was first introduced by Wiener in 1924 [2^, in which he showed that the probability distribution 
for a classical system is found by summing the final states of all paths that could describe the evolution of the system. 
The corresponding result for quantum systems was proven by Feynman in 1948 j29j , and a version that uses the 
Hamiltonian, rather the Lagrangian, is given by Farhi and Guttman [30| . The latter approach is used by Mackenzie, 
Marcotte and Paquette (MMP) [2^ to find an expression for the time-evolution of adiabatic quantum systems using 
perturbation theory. We provide the background material that is needed to understand what we mean by paths in this 
section, and in the subsequent section, we use these results to derive the MMP representation of the time-evolution 
operator, but without using perturbation theory. 

Our usage of path integrals is somewhat different than Feynman's original formulation. We obtain our path 
integral representation, similarly to Feynman, by first discretizing time. By approximating the Hamiltonian with 
one that is piecewise constant on these discrete time steps, we obtain an approximate evolution operator that is a 
product of exponentials of the Hamiltonian evaluated at each of the discrete times. Unlike Feynman, we switch to 
the instantaneous eigenbasis of the Hamiltonian at the given time-step. The resulting expression is then simplified 
by expanding the product of sums into a sum of products of the projectors that are used to represent the basis 
transformations. Each product of projectors that occurs in this expansion can then be interpreted as a description 
of a potential evolution of the system, which we call a path. We then take the limit as the number of time-steps in 
our discrete time-evolution operator approaches infinity. This removes the approximation error caused by taking the 
Hamiltonian to be piecewise constant and turns the sum over all such paths into a sum of path integrals. These ideas 
will be elaborated on as we proceed. 

We now present a more formal discussion what we mean by a path. A path is a description of a sequence of states 
that could describe a sequence of observations of the evolution of the quantum state. In particular, we take a path 
to be an ordered pair {{i'}, {s}), where {ly} is a list of labels that describe the instantaneous eigenstates of II{s) and 
{s} specifies the time that the system transitions to that state. 

We require that no two consecutive labels in a path are identical, ie., Vq ^ i^q+i for all q. With this restriction, it 
is clear that the elements of {s} correspond to the times at which the system jumps from one state to another. We 
also take the sequence s to be a strictly increasing sequence of dimensionless times taken from the interval [0, 1]. We 
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FIG. 1: These diagrams are a representation of the paths used in the examples in this section. Here the state labels refer to 
different values of v and do not necessarily refer to energy. 



take these times to be strictly increasing because paths that move backwards in time, or are in two different states 
at the same time, cannot describe a physically realistic sequence of measurements of the state in its instantaneous 
eigenbasis. As a final note, we first consider the case where each path is composed of at most L+1 distinct states and 
times, where the times in {s} are taken to be integer multiples of 1/L to simplify the analysis, similarly to Ambainis 
and Regev [25j . Unlike the formulation of their result, we take the limit as L approaches infinity early in our analysis 
thereby removing L from our final results. 

As an example of a path, consider the evolution for a three-state system with L = 4 given in Fig. [T^: |0(0)) — > 
|l(l/2)) — ^ This sequence of states can be recorded as a path that is initially in the state |0(0)) and jumps 

to the state |l(l/2)) at time s = 1/2. The resulting path ({i^},{s}) is given by the sequences {v} = {0,1} and 
{s} ~ {0, 1/2}. Equivalently, this path can be expressed as a diagram, as seen in Fig. [T^. 

Diagrams are a convenient way to think about the paths that describe adiabatic evolution. In our diagrams, an 
arrow denotes the evolution of the state as a function of time. If the arrow runs along a line, which denotes a particular 
energy level, then it implies that the system remains in that instantaneous eigenstate throughout that portion of the 
evolution. If the arrow connects two energy levels, then it implies that a transition, or jump, between energy levels 
occurs at that time. As level crossings are forbidden in our analysis, we do not change the relative spacing between 
the energy levels as the eigenvalue gap changes with time. Such diagrams are useful because they convey the intuition 
of what is being described by our more mathematical notation, which uses a tuple of sequences to describe the path. 
A similar diagrammatic notation is also presented by MMP [26]. 

In this work we classify paths by the number of jumps that they contain. As an example, consider the following 
evolution: |0(0)) |l(l/2)) ^ |2(3/4)) |0(1)). The corresponding path is parameterized by {v} = {0,1,2,0} and 
{s} = {0, 1/2, 3/4, 1}. This path jumps from one level to another three times and therefore we say that it contains 
three jumps. Similarly, the previous example contained one jump. The number of jumps present in a path is \{v}\ — 1. 
If {{v}, {s}) is a q-jump path then we define {v} :— {i/q, vi, . . . , Vq\ and similarly {s} := {sq, si, . . . , Sq\, where vq = G 
and sq = 0. We choose vq and sq in this fashion so that the p"^ jump occurs at time Sp, which makes our results 
easier to read. A diagram of the path in this example is given in Fig. [T}d. 

Additionally, some paths describe adiabatic evolution and others do not. We call an n-jump path "adiabatic" if 
the final and initial state in the path are the same. For example, the path in Fig. [Jd, is an adiabatic path whereas 
the path in Fig.[T^ is not. We call any path that is not "adiabatic" a "non-adiabatic" path. We use the notation Jq 
to represent the set of all g-jump paths and to represent the set of all (/-jump non-adiabatic paths. Non-adiabatic 
paths are crucially important to our work because the error in the adiabatic approximation is given by the sum over 
all non-adiabaticpaths. This observation is especially relevant in light of the fact that the convergence problems 
noted by MMP [2^ occur because their sums are not restricted to non-adiabatic paths. 

We describe below, the sum over all such paths in the limit of continuous time. 

Definition 2. Let ({i'},{s}) represent an arbitrary path in Jq, in the limit as L oo. We then define the sum over 
the set all q-jump paths, Jq to be 

dW:=E E--- E / / (11) 

/ I J J Jo 
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and similarly we define the sum over all q-jump non-adiabatic paths, J'^, to be 

E E E / •••/ / (12) 

where ds* — dsi • • • dsq. 

Now that we have introduced our notions of paths and the sum over all paths, we apply them to describe the 
time-evolution of quantum mechanical systems. 

V. PATH INTEGRAL REPRESENTATION OF TIME-EVOLUTION 

In this section we show that by taking the continuous limit of a discrete approach that is similar to that of Ambainis 
and Regev [25| . we obtain a path integral representation of the time-evolution operator that is similar to that of 
Mackenzie, Marcotte and Paquette (MMP) [2g|, but without using perturbation theory. The convergence issues noted 
in [2^ are ameliorated by noting that not all n-jump paths contribute to the error in the adiabatic approximation; 
however, the sum of all n-jump non-adiabatic paths does. Our goal in this section is to find an expression for the 
sum of all the non-adiabatic paths to the error in the adiabatic approximation. Upper and lower bounds for the 
contribution of these paths are proven in Sec. I VIII 

To simplify our proof, we choose the eigenstates of the Hamiltonian such that 




I if En{s) = E,n{s) 

Pn,m{s) := {h{s)\m{s)) = { (n(s) \H{s)\m(s)) otherwise ^^^^ 



similarly to [25l [26j . Here Pn,m{s) can be interpreted as a tunneling amplitude per unit of dimensionless time 
between the states |n) and |m). The choice /3„ „(s) = corresponds to choosing the phases of the instantaneous 
eigenstates (which are dynamically irrelevant [sij ) such that the geometric phase [Tsj is zero, which simplifies our 
proof dramatically. We use this choice of eigenstates to state in the following theorem, which is the main result of 
this section. 

Theorem 3. //||_ff(s)|| is bounded and H{s) is piecewise- differentiate on the interval [0,1] then ?7(1, 0)|G(0)) is 
C/(1,0)|G(0)) = 

e-/o^i^o(Od«T|g(^)^^ j (e-<^-«)^«+/o'^-(«)'^«)^/3,,,G(si)) |^^i(l))d(Ji) 
+ j (e-^(/4^-(W=?^-(Wo^^-(«)<^«^^^^^ |.2(l))d(J2) + --- , (14) 

where {{v}t{s}) represents an arbitrary path in Jq and Vi and Si refer to particular element in the sequences, and (3 
is defined in Eq. 

The series in Eq. [TJ] has an intuitive interpretation. It states that the time-evolution is a weighted sum over all 
paths that could describe the time-evolution of the quantum system, wherein the weight of a particular path is given 
by the products of the tunneling amplitudes multiplied by the dynamic phase picked up in that path. Adiabatic 
behavior occurs because in the limit of large T, the phases that are achieved by each path connecting |G(0)) and 
|i^(l)) for v G, vary considerably. Similarly, we expect that the contribution of the sum of all paths that begin in 
|G(0)) and end in |G(1)) to be suppressed given that the path contains two or more jumps. Hence we expect that 
interference effects will cause their contribution to become diminished in this limit, leaving only the zero jump path 
which links |G(0)) to |G(1)). This phase cancelation effect is demonstrated in Fig. [5] 

We use these properties of the eigenstates to show the following lemma, which is used to evaluate the products of 
projectors that appear in our proof of Theorem [31 

Lemma 1. Let {{i'}, {s}) be a path in Jq with {s} — {0,ji/L, . . . ,jq/L} for an integer value of L then, 





Mq/L)){vm (^nS^^±i^ + o(jVi^+^)) , (15) 

where each P^.(j/L) = / L)) / L)\ projects onto the ray spanned by the state \vi{j / L)) . 
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FIG. 2: This figure is a plot of tiie phase of states output by 21 one-jump paths that jump at time si = {0, 1/20, 2/20, • . • , 1} 
from the ground state and the first excited state for a Hamiltonian with Eg{s) — and -Ei(s) = 1 + s, for T — 0.01, 
T — 1 and T = 4 in a, b and c respectively. Here each solid line is a unit vector in C that is oriented in the direction 
of exp(-i X, Ei{^) - EG{0<i^T). This plot shows that as T increases the vectors tend to become uniformly distributed 
throughout the unit circle, suggesting that the sum over all paths becomes less substantial as T increases. 



Proof. We prove the lemma by performing induction over q, which we use to denote the number of jumps in the paths 
in question. Our first step is to prove the base case, q — 1. Using Taylor's theorem, and the fact that the eigenstates 
obey (lyilvi) — 0, it follows that for any i^i and s, 

(Ms + 1/L)\M.'^)) = 1 + 0{l/L^). (16) 

This implies that for any 

f[ P.Aj/L) - \M3i/L))Wi{[ji-i + 1]/L)\ (1 + O ([j, - - 1]/L')) . (17) 
i=j<!-i+i 

We demonstrate the base case, g = 1, of Eq. [T5]by substituting £ = 1 and jo = into Eq. [T71 
Now let us assume that the result of Lemma [1] applies for the case 



n ^«.-.fi;) - n« 



U=J,-2+l ^ ^ / Vj = i ^ 



= |.,_i(,,_i/L))(.i(0)| (lMdj^l±p^ + o{j,-,/L'')^ . (18) 
By using Eq. [TS] and Eq. [T7] we find that the left hand side of Eq. [12] becomes 

Taylor's Theorem and the fact that, from the definition of a path, Vg =/= Vq-i imply, 

+ l]/L)|z.,_i(j,_i/i)) = i(;>,(j,_i/i)|z.,_i(j,_i/i)) + 0(1/^2). (20) 

^ /3.,,.,-,(j,-iA) ^^^^^^2)^ (21) 
L 

We substitute Eq. [5T]into Eq. [T21 The induction step on q follows, completing our proof of the lemma. □ 
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Lemma [T] provides an expression that can be used to evaluate products of projection operators. These results allow 
us to prove Theorem [3] by discretizing the time-evolution and examining the evolution in each time step in the basis 
formed by the instantaneous eigenstates of H{s). We use projectors to transform the basis from the previous time 
step to the basis at the current time step. The resulting expression for the time-evolution is given as a weighted sum 
of products of projectors, which we then simplify by using Lemma [T] Our proof of Theorem [3] is given below. 

Proof of Theorem\^ It follows from using the Trotter formula [s^ for a fixed value of L that if H{s) is piecewise 
differentiable then we can approximate H{s) by a piecewise constant Hamiltonian to find 



L /N-1 

= n E 

j=0 \iy=0 



^^•^(J/L^T/L + 0(1/L), (22) 



where for any s € [0, 1] and v € {0,1, . . . , N — 1}, \v{s)) is an eigenstate of H{s) and Pp{s) := \i'{s)){i'{s)\. We then 
expand Eq. [21] to find 

C/(l, 0)|G(0)) = E • • • E n ^'"-^^"-P., 07i)|G(0)) + 0(1/L). (23) 

Equation [531 is a weighted sum over every possible path that can describe the evolution of the quantum system over 
the time interval. We group the paths in the sum by the number of jumps that comprise them and apply the result 
of Lemma [1] to find, 



C/(l, 0)|G(0)) =6-^(5:^=1 iJo07i)T/L|g(^)^ 



L 

e 



E E z ^--^(^^/^V i(i)) 

+ E E ^ ^^^e~'^^^=^2 B„,07L)+E^:L-'£;Gb7i 

vii^vx ui^G j2=2 ji = l 

l3u2.vl{j2lL)Pi,^,G{illL) \ 
X J2 \^2[1)) 

+ --- + 0{l/L). (24) 

Equation [Ml is similar to that of Theorem [31 except instead of time- integrals we have discrete sums over time and we 
also have an 0{1/L) term that represents the error in simulating the adiabatic evolution using a piecewise constant 
Hamiltonian. We replace many of the sums in Eq. [24lby integrals, using the fact that 



J2 fij/L)/L = r f{s)ds + 0{l/L), 



(25) 



j=siL 



and eliminate the simulation error, by examining Eq. [24lin the limit as L ^ cxi. This demonstrates that the sum of 
all one- and two-jump paths approach the path integrals in Theorem [31 The analogous results for paths that contain 
three or more jumps follow from the exact same reasoning, completing our proof of the theorem. □ 

We have introduced our concept of path-integrals in this section, and have found an expression for the time-evolution 
operator in terms of these path integrals. Our next step is to approximate the contribution of the non-adiabatic paths 
in Eq. 1141 This gives us the error in the adiabatic approximation. We use the results of the following section to find 
such a bound, which we give in Sec. I VIII 



VI. A RECURSIVE EXPRESSION FOR PATH INTEGRALS 



Although an upper bound on the eontribution one-jump non-adiabatic paths to the error in the adiabatic approxi- 
mation can be directly calculated from Theorem[31 it is difficult to upper bound the contribution of paths that contain 
more than one jump using that strategy. To do so, we use the fact that every m-jump path is a concatenation of 
m one-jump paths. We use this property to find a recursion relation that allows us to express the contribution of 
an m-jump non-adiabatic path in terms of the contribution of an (to — l)-jump non-adiabatic path (this fact follows 
from our definition of paths that we introduce in Sec. IIV|) . Our error bounds then follow from summing a geometric 
series that results from unfolding these recursion relations. 

Using the series in Eq. [14] and neglecting all paths that end in G (because such paths do not contribute to the 
approximation error), the error in the adiabatic approximation can be expressed as 

((]I-PG(l))C/(l,0))|G(0))=f; j ^ (1)1^.(1)) ) ' (26) 

where 7^,1/(5) represents the eigenvalue gap between the instantaneous eigenstates |/i(s)) and \v{s)) and /jy^ represents 
the sum over all times as given by 

X / 'e-*^o^''o..,(?)dCT^^^_g(^^)dsi...ds,. (27) 
Jo 

The integrals in Eq. [57] are expressed recursively as 

US^) = re-^^o^^''--'(«''^«^/3.,,.,_,(s)/.,_,(s)ds, 
Jo 

UAx)^ re-'/o^'^-i(«)'^«X,G(s)ds. (28) 

JQ 

A similar recursive expansion to Eq. [55] is found by Sun in 27]. In the following section we find higher-order 
approximations to the error in the adiabatic approximation by using integration by parts to approximate the recursion 
relations in Eq. [28] 



VII. APPROXIMATING PATH INTEGRALS 



We have argued qualitatively that quantum interference effects are responsible for the emergence of adiabatic 
behavior in quantum systems. In this section we provide quantitative bounds for this phase cancelation effect. We 
first provide an example of how integration by parts can be used to approximate the phase cancelation effect for 
one-jump paths, and then use the same strategy to build an asymptotic expansion for the sum over all g-jump non- 
adiabatic paths in powers of T~^, for every g > 0. We then sum these expressions over all values of q and obtain 
bounds for the approximation error. We then prove Theorem [5] using these bounds for the error. 



A. Approximating Integrals Using Integration By Parts 

In this subsection we provide an example that shows how integration by parts can be used to estimate the contribu- 
tion of all one-jump non-adiabatic paths to the error in the adiabatic approximation. The purpose of this introduction 
is to provide a simple example of approximating path integrals that will help explain our approach for bounding the 
contribution of q-jump paths to the approximation error, which we address in the following subsection. 

Integration by parts is a technique that is commonly used to approximate the integrals that appear in proofs of the 
adiabatic approximation [l^, H^. Unlike other methods, such as the method of steepest descent, integration by parts 
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gives an exact expression for the integral. Specifically, if / : 
dsg{s) is non-zero in the region of integration then [33| 



i-^ C and g : R I— >■ C are differentiable functions and if 



ds 




ds 



(29) 



Equation [52] gives the leading order behavior of the integral in the limit of large A as well as an exact expression 
for the remainder term (which can also be approximated using integration by parts) given that dsg{s) ^ for all 
s G [0,x]. For this reason, we assume that the minimum eigenvalue gap between any two non-degenerate eigenvectors 
is non-zero because, if it is ever zero, then integration by parts will fail to approximate Eg. 1141 In such circumstances, 
an alternative method would have to be used to analyze the integral. 

We approximate the contribution of all one-jump paths to the error in the adiabatic approximation as an example 
using integration by parts to bound a path integral. The contribution of all one-jump paths to the approximation 
error is 



|i^i(l))/3.,.G(s)e-*/o^o,.,(«)d?T^(j/) 



: ^ |i.i(l)) / /3,„G(s)e-^^o^c=,.i(?)d«Tj^_ 



(30) 



We apply integration by parts to obtain, 



1^1 #G 



s=0 



d_ 

ds 



-ilG,uAs)T 



(31) 



Using :— niaxs ||i7(s)|| and maxs ||-ff(s)|| G 0{l) we find that Eq. [3T]has the upper bound 



\v^{l))(i,,^G{s)e 



-*/o°7G,, 



< 2- 



\H\\ 



sb.^,G{s)?T 



0[l/T% 



(32) 



which is a commonly used estimate of the error in the adiabatic approximation. If the contribution of all non-adiabatic 
paths with two or more jumps is negligible, then Eq. |32] implies that the standard estimate for the approximation 
error holds [26*1. 

Higher-order estimates for the contribution of the one-jump terms can be found by applying integration by parts 
to approximate the remainder term. The success of this strategy depends on the radius of convergence of the 
perturbation expansion. Unfortunately, the series has a zero-radius of convergence. This implies that arbitrarily 
high-order adiabatic theorems generated using our technique will only be useful in the limit as T — oo. This is 
because the repeated derivatives that occur from using integration by parts to obtain an 0{1/T*') approximation to 
Ci are 0{k\/T^). The perturbation series for the one-jump path integrals will generally have a radius of convergence 
of zero as fc — >■ oo; suggesting that higher-order adiabatic approximations are often only useful in the limit of large T, 
as also noted in (l3| . 

In the subsequent subsection we generalize the ideas presented in this subsection to find bounds for the approxima- 
tion error that emerges due to every non-adiabatic path, as opposed to the results in this subsection, which considered 
only paths with one jump. 



B. Error Bounds For The Adiabatic Approximation 

Our strategy for bounding the error in the adiabatic approximation is to use integration by parts to find upper 
bounds for the oscillatory integrals present in Theorem [3] Our zeroth-order and first-order adiabatic approximations 
in Theorems [1] and [21 directly follow from this upper bound, which is given in the following Theorem. 

Theorem 4. Let H : [0, 1] i— > C^^^ he a Hamiltonian that is differentiable three times and has a minimum eigenvalue 
gap between any two non-degenerate eigenvectors of 7min > 0. Then, for any T > 0, the error in the adiabatic 
approximation obeys 



(i?-PG(l))C/(l,0)|G(0))- E 

i^^G 



e-*^"|z/(l))(z>(si)|G(si))e-*/o'^G..(«)d«T 



-ilGAsi)T 



si=0 



<- 



R 



(33) 



11 



where 4>i, = /J^ i?^(^)d^r and R is hounded above by 

.211^11 + llijll 25||ij||||H|| + 16||ij||2 



R<- 



\HP 12||if||||ii-||2 + 118||iJ||3 36||ij||4 



73. 

' min 



1^- 

I mm 



75. 

' mm 



76. 

/ mm 



r 



and V — 



\\H\m\ 

— 



\\H\\T 

2\\H\\^ 



r/(7„ 



.T) _ 1) _ 



T 



(34) 



-ff II , ||i?|| and \\H\\ refer to the maximum values of the norms of the first three 



Tmin / Tmic 

derivatives of H and 70,1^(5) = Eg{s) — E,j{s). 

This theorem is important because it provides upper and lower bounds that, except in pathological cases, converge 
to the observed approximation error in the limit of large T. As a result, these bounds are superior to the JRS 
bound [13] in the limit of large T. We compare these two bounds in greater detail in Sec. I Villi 

The proof of Theorem 2] requires several steps, which we express in three lemmas. In order to abbreviate the 
statements of our Lemmas, we define the contribution of all j-jump non-adiabatic paths to the error in the adiabatic 
approximation to be Cj . Hence the error is 



where for example. 



(/ - PG(l))t/(l, 0)|G(0)) = Ci + C2 + C3 + 



(35) 



(36) 



and 



Co 



(37) 



We have that 1^2 7^ G in our expression for C2 because any such path is an adiabatic path and thus does not 
contribute to the error. The following lemma contains our upper bounds for the contribution of Cj for j > 2, which 
is an important step towards proving Theorem 

Lemma 2. Let H : [0, 1] 1— > C^^^ be a Hamiltonian that is twice differentiable, then the following bound holds on 
the error in the adiabatic approximation: 



\{1 - Pg(1))C/(1,0)|G(0)) - Cill <||C2|| + 1 + 



r 



|i?llT 



r/(7mi„T) 



(38) 



Here F 



\\H\\\\H\ 



mi 

1, 



■j-^ and \\H\\ and \\H\\ refer to the maximum values of the norms of the first two 
derivatives of H and jmin is the minimum eigenvalue gap between any two non- degenerate eigenvalues. 

Our proof of Lemma [5] uses integration by parts to approximate ||Gj|| for j > 2. To do so we must provide upper 
bounds for the derivatives of /? that appear due to using integration by parts to evaluate the integrals. The most 
straightforward way to approximate these derivatives is to differentiate Eq. [T2] and use the triangle inequality on the 
result. However, this simple approach leads to extraneous factors of N, which appear because of the sums over i/q. 
We use a trick involving the resolution of unity to provide upper bounds for the first two derivatives of /3 that do not 
depend on N. We state this result in the subsequent lemma, and will use it to prove Lemma [H 

Lemma 3. Let H{s) be twice differentiable on s € I = [0, 1] and let s' G /, x : Jg_i C &e a bounded function, and 



for simplicity we define the argument of x to be the path v :— {{vj : j — 0, . . . ,q ~ 1}, {sj : j = 0, 



1}) the 



and 



E /3^,.^,-i(^)x(^)k,(5')) 



l/„M„^l,...,iyi 



< 



i/<,,i',_i,...,i'i 



< 



73. 

' mm 



4||i?||^ . IliJil 



2 

7min 7min 



E x(^)i^,-i(s)) 



I^„_l,l/„_2,...,l'l 



UWHf , 12||iJ||||ffl| , IlijII 



72. 

/mm 



7n 



E X(^)l^,-l(^)> 



i/g-l ,I/<j_2,...,I'l 



(39) 



(40) 
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The function x in these bounds is an arbitrary function, for example if we wish to find an upper bound for 
Sjyj vi Pv'i,ui{^)Pi'i,G{s)\v2(s')) , then we can use Eq. |40]to find the upper bound by choosing x(i^i,s) — I3^^,c{s)- 
The proof of Lemma [3] is given in Appendix \X\ 

We will now use the first result in Lemma [3] to prove Lemma [21 

Proof of Lemma [H Using integration by parts on Eq. [2H1 yields 



d ( P>^^,y,-i{s)fy^-As) 



-.^,-i.^,(-)rds, (41) 



where Ku^_^^u^{s) = 7i/,_i,iy, (C)dC and fug{x) represents the time-integrals in the sum over all g-jump non-adiabatic 
paths. Using Eq . and the Fundamental Theorem of Calculus, we can simplify Eq. UTJby substituting 

dsf.g^A^) = e-'''-.-''.-(^'^/?.,_,^.,_,(.)/.,_,(.). (42) 

We also have that for any p > 0, /p(0) = 0. We then simplify Eq. |4T]by using the recursion relations to express every 
instance of as an integral of fvc-2^ to find 



^.,-.,.,-1 (^i)^;3^^_^^^^_^ (si)/.,_, (Sl)dsi 



— / e-'='^.-'^.-(^^)X-i.^,-.(^2)/.,_.(s2)d2 



9*1 \-ilu,-i,u,{si)T , 

1 V -«7i',-i,i^,(s)r / 
We then express the sum over all values of z^i , . . . , z/g in Eq. 23] as 



(43) 



(44) 



where for any q > 0, 



Vq,Vq-l,... 



f q ,f q — 1 ,■ ■ ■ 



9si -i7,.,_i,i.,(si)r^ 

X /3,.,_i,i.,_2(s2)/t.,_2(s2)d2s, 



(45) 



where d^s = ds2dsi and 



-«7i',-i,i',(a^)r 

(^1)/^,-. (^1) 



-j7>^,_i,i',(a;)r 



dsi 



(46) 



We now will proceed to bound ||(7g(a;)||. To do so, we use the definition of /3 in Eq. [T3]to find following bound, 
which is valid for any s' , s" , s'" and s"" . 

E l^,(^))/3.„.,_.(^^')/3.,-.,.,-.(^")/.,^.(^"')/7.,-...,(s"") 

E.„.,_„...l'^,(^))(^,(«')|i?(^')l^.-l(s'))K-l(^")|ij(s")l^.-2(s"))/.,_.(s"') 



"1Vq-l,Vq (•S')7l'g-2,l'g-l i^" )^l^q^ 1 ,l^q (s"") 



< 



< 



73. 

' mm 



E ^,-.(^"')l'^,-2(s")) 



E /'',-2(s"') 1^^9-2(5")) 



(47) 
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where 7min is the minimum eigenvalue gap between any two non-degenerate eigenstates, for any s €E [0, 1]. We then 
have 



'111 in 

2\\HP 



E 



dsi. 



< 



73. T 

I mill 



(||5g-2(si)|| + ||V2(si)||)dsi, 



(48) 



where the last inequality follows from the triangle inequality. Using similar reasoning and our expression for the 
derivative of /3 in Eq. [39l we find that 




^,_2(S2)|J^5-2(S2)) 



(||ff,_2(si)|| + ||V2(^i)ll)d^s 



(49) 



We eliminate |lgg-2|l and ||/ig-2|| in Eg. H51 and Eq. iterativelv by using Eg. H51 and Eq. to provide upper 
bounds for |j(7g_2|| and ||/ig_2|| in terms of ||5q-4|| and ||ft.q_4||. This process is then repeated [<7/2j times to eliminate 
\\hi\\ and \\gi\\ from our bounds for any i. In the case where q is even, no integrals will remain after performing 
this substitution q/2 times, whereas if q is odd, the remaining terms will have one integral left that has not been 
approximated using integration by parts. 

Let us now consider the case where q = 2rn. We expand fu2mJ which is the time-integral part of 6*2™, into a sum 
of 2™ terms by applying the recursion relations in Eq. |49] and Eq. 1481 We then group these terms by the number of 
integrals present in a given term. It follows from combinatorics that, for any i e {0, 1, . . . , to}, there are (™) ways to 
obtain a term with 2to — £ integrals and hence 



\\h2rn{x)\\ + \\g2m{x)\\ 



< 

i=0 

< 



2HV\ H H 6 Lff 



1 4^ rm\ f 2\\H\\A"'''^ f \\H\\\\H\\ 6\\HP^ 



■E 



TO!(7mi„T)-^ V ^-i" / V ^-i" ^= 



TO!(7„i„T)™' 

where 1/(2to — i)l is the volume of integration and the last line of Eq . [5(11 follows from the binomial theorem for 

^. f2\\Hr + \\H\\\\H\\ , 6||i7f \ 



^min ^iiiin 



(50) 



(51) 



By summing Eq. 1501 over all to > 2, we find that 



Er2„.|l<er/(>-"^)--i--l. (52) 

o 7niin-^ 
m— 2 

Our next step is to find J2m=i l|C'2m+i||- To do so, we use the fact that a path that contains 2to -|- 1 jumps is 
the concatenation of a path that contains 2to jumps and a path that contains one jump. An upper bound on the 
contribution of non-adiabatic paths with 2to + 1 jumps is found by taking the norm of the result for the 2to jumps 
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in Eq. [50] and multiplying it by an upper bound on the remaining one-jump path. Hence an upper bound on the 
contribution of these paths is 



\C: 



2m+l\ 



< 



< 



m!(7„,i„T)™ sG[04] 

pm 



< 



d_ 



-ilvi,G{s')T ^ 

r 



-ilutfi{s')T 



m!(7mi„T)™ \\H\\T 
By summing this result over m = 1, 2, . . . we find that 



< 



rn — l 



\H\\T 



„r/(7„inT) 



Adding Eq.[52]to Eq. [Ml yields 



Eiic.ii< (i 



J=3 



,r/(7„.„T) 



r 



7minT 



(53) 



(54) 



(55) 



The proof of the Theorem then follows from the above results. By definition, the contribution of the one-jump and 
two-jump paths to the error in the adiabatic approximation are ||Ci|| and IIC2II respectively. Equation [351 then follows 
from the triangle inequality. □ 

Now all that remains in our proof of the adiabatic approximation is to bound the one- and two-jump contributions 
to the error in the adiabatic approximation. The following lemma contains our bounds for ||Ci|| and IIC2II. 



Lemma 4. Let H 



NxN 



be differ entiahle three times then, 



-iK.G,u{si)T 



where 4>i, = Ei,{^)d^T, 



Rn 



2\\H\\ 



H 



73. T2 

I mm 



20||i?||||iJ|| + 12||i7|| 
74. T2 

/ mm 



<Ro, 



(56) 



\H\ 



rp2 ' 



(57) 



IC2II < 



\H\\ 



?>\\H\\\\H\\ 



mmmw 



30||i/|| 



74. r2 

/mm 



7^. T2 

/ mm 



36pr 

76. T2 

/ mm 



(58) 



where \\H\\ — max^g/ (s)|| , \\H\\ = max^g/ (s)|| , = max^g/ ||_ff(s)|| anrf 7inin is the minimum energy gap 
between any two non- degenerate energy levels. 

The proof of Lemma [4| is similar to our proof of Lemma [2] The main difference between our proofs is that Eq. [43] 
is not used, and instead we apply integration by parts twice to guarantee that our approximations are correct to 
0(l/r2). We then use the triangle inequality to find upper bounds for the 0(1/T2) terms and then use the triangle 
inequality to prove our error bounds. The details of this proof are given in Appendix [Bl 

The proof of Theorem [4| follows from the results of Lemmas [2| and [4] 

Proof of Theorem^ The proof immediately follows by substituting the bounds for ||Ci|| and IIC2II in Lemma [4| into 
the result of Lemma [2] □ 
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The proof of our zeroth-order adiabatic theorem, given as Theorem [Tl follows directly from this result. Similarly to 
the proof of Theorem 21 the proof of Theorem [5] follows from combining Lemma [5] and Lemma S] and then showing 
that the terms that are proportional to T^^ are 0(Ai/7i„inr^), where the timescale Ai is given in Definition [T] 

Proof of Theorem\^ Our asymptotic bounds are taken in the limit of small Ai/[7niinr^] and large Ai. From the 
assumption that Ai is asymptotically greater than any constant, we find by substitution that each of the five terms 
that are in Eq. [56]are proportional to and their sum is 0(Ai/[7minr^]). This implies that 



1 

+ 0(Ai/[7„,inT2]). (59) 







-ilGAsi)T 

It follows from substitution into Eq. [58]that C2 G 0(Ai/[7niinT'^]). 
Similarly, we have that 

\ ^min 7min 7min / 

The second term in the above equation is asymptotically greater than or equal to the first term if || G 0(||ij|p/7min)- 
However, this condition implies that the third term in Eq. [60]is asymptotically dominant. Therefore, 

Y '^min ^min / 

The first term in the above equation is dominant only if ||i^||^/7min ^ It follows from Eg . [6T1 that 

r2 e o j + j e 0(7minAi). (62) 

From our assumption that Ai/(7inin2^^) ^ wc find, 

e (63) 

/min 

which implies that 

G 0(1). (64) 

This result is important because it allows us to replace exp(r/[7niinT]) with a truncated Taylor series, without 
introducing any error terms that are exponentially large. Following this strategy, we find from using Taylor's the- 
orem and Eq. [64] that the result in Eq. [55] is 0(Ai/[7ininT^]), proving that the contribution of HCsH, ||C4||, . . . is 
0(Ai/[7ininT^]), which proves the theorem. □ 

We then can state the proof of our zeroth-order adiabatic theorem as follows. 

Proof of Theorem [IJ The proof directly follows by applying the triangle inequality to the result of Theorem [5] and 
maximizing each term that is present. □ 

Our upper bounds can be further tightened if we assume that the Hamiltonian has only two distinct energy levels 
for all s. This is because, in this case, all paths that contain an even number of jumps and start in the state |G(0)), 
must end in the state |G(1)) because of our choice of eigenvectors in Eq. [13] Therefore we neglect all paths that 
contain an even number of jumps. This result is stated in the following corollary. 

Corollary 3. Let H{s) he a twice dijjerentiable Hamiltonian with two distinct energy levels Eg{s) and Ei{s), where 
the eigenvectors corresponding to Ei{s) are (N^l)-fold degenerate. Given these assumptions, the error in the adiabatic 
approximation is bounded above and below by, 

Wil - Pg(1))C/(1,0)|G(0)) - Cill < (^-J— ) (e//(^™^) - 1), (65) 
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Proof. The proof of this corollary follows from our proof of the previous theorem. Although Ei is {N — l)-fold 
degenerate, this system is an effective two-level system because our choice of eigenstates prevents any transition 
between two states with the same energy eigenvalue. Given this choice of eigenbasis, any path that has an even 
number of jumps is an adiabatic path, which are excluded from the sum. 

Equation [54] then gives an upper bound on contribution to the error due to paths with more than one-jump because 
all paths with an even number of jumps do not contribute to the error in this case. Equation [65] follows from the 
definition of Ci and the triangle inequality. □ 



VIII. SEARCH HAMILTONIAN 



In this section we _apply our error bounds in Theorem[3]to the case of a Search Hamiltonian, and compare our results 
to the JRS bound [10| . In addition, we discuss a new method for choosing the integration time steps geometrically. 

The Search Hamiltonian is a Hamiltonian whose ground state at s = is an initial guess for the solution to a 
search problem and the ground state at s = 1 is chosen to be the state that is sought by the algorithm 0, [HI . If the 
Hamiltonian evolution obeys the adiabatic approximation then the initial guessed state will be transformed into the 
final solution state, thereby solving the search problem. 

We denote the initial state jV'o) and the marked state \ipi), and the goal of the algorithm is to choose a Hamiltonian 
that generates an evolution that transforms the initial guessed state into the marked state in the adiabatic limit. We 
evolve our system according to a time-dependent Hamiltonian that is a linear interpolation between these two: 

H{.3) = (1 - s)Ho + = a - (1 - s)|V'o)(^o| - sIV'iX^il- 

Although the search Hamiltonian has N eigenstates, only |l(s')) has non-zero overlap with the ground state \G{s)) 
for any s and s' . Therefore the search Hamiltonian generates an effective two-dimensional evolution. The eigenvalue 
gap between for the two relevant states, if l-^o) = '^x=o |2^)/V^ IV'i) is a state in the computational basis, is [ll| 



7i,G(s) = yi-4(^l-l^s(l-s). (66) 
If we define 9 — arccos(|('0o|''/'i)l) then the instantaneous eigenvectors of the Hamiltonian are 



and 



where 



|l(^)> = ^^^^^^^l^o)-^|^r), (68) 
sm V sm 



. N 1 f ssin(26l) , 
ip{s) = - arctan ^ — . (69) 



2 V(l - s + scos(26l)) 

These expressions for the eigenvectors clearly show that the eigenvectors undergo rotation in the plane formed by 
the span of jV'o) and \ipi), where the rotation angle is given by (p{s). Similarly, the rotation speed is given by 93(5). 
This implies that the eigenstates rotate more slowly near s = 1/2 than near s = or s = 1, it makes sense to use an 
integration scheme that takes constant sized steps in ip rather than constant sized steps in s. This is because constant 
sized steps in s will lead to overly conservative time steps being taken near s = or s = 1. 
We find for the case where 9 = arccos(l/\/]V) that 

|"arctan(VA^ - l(2s - 1)) + arctan(ViV - 1)1 , , 

<p{s) = ^ ^ ^ ^. (70) 

Interestingly, this parametrization has been used in the literature to interpolate adiabatic evolution in order to 
minimize the approximation error (ill . |33 |. To clarify, we use this interpolation for the timesteps taken in our 
numerical method, but use linear interpolation for the Hamiltonian. 

We use this result to test our upper and lower bounds for the error in the adiabatic approximation by computing the 
time-evolution numerically using a modified version of Eq. [21] that uses evenly spaced steps in if rather than s. The 



17 



10-' 



Q. 

E 

< 10-' 
i_ 

2 

UJ 



10-' 



10-" 



^□□□□□BDD 



— s — Lowest Order 
□ Our Upper Bound 

JRS Bound 

Numerical Result 




10" 



10' 



10' 



10^ 



10" 



10" 



10° 



10' 



FIG. 3: This figure sliows tlie excitation amplitude for a search Hamiltonian with 4 potential marked states with Itpo) 
X]J!=o \^)/^^- The jaggedness of the figure appears because of rapid oscillations that occur in the limit of large T. 



results were found by adaptively choosing L such that doubling its value resulted in less than a one percent change in 
the computed error in the adiabatic approximation. To test the results, we also compared the results of the simulation 
for small values of T to those found by solving the Schrodinger equation using an adaptive quadrature Runge-Kutta 
integrator using double precision data and found excellent agreement. Ho-wever, we found that the Runge-Kutta 
integrator was unable to provide sufficient accuracy for simulations with large T and hence we used Eq. [52] in its 
place. 

In Fig. [3] we use these expressions for the eigenvalue gap as well as the eigenvectors to evaluate our upper bound 
in Corollary [Sj and compare them against the JRS bound [lol | as well as our numerical data and our lowest order 
approximation to the error as given by Theorem [5] In all cases, our upper bound on the error is superior to the JRS 
bound, in the limit of large T. We also note that the JRS bound does not converge to the observed error in the limit 
of large T, whereas ours does. As expected from the fact that our bounds are higher-order than the JRS bound, our 
bounds are inferior for small values of T. 

Fig. m demonstrates that our upper bound and lower bounds approach the result of Corollary [3] for large T. We 
do not plot our numerical data here because for T > 100 the numerical results are graphically indistinguishable from 
the result of Corollary [3] and the computation of the error in the approximation becomes prohibitively expensive for 
T > 10^ 

We cannot expect our bounds to be asymptotically tight for the search Hamiltonian because the leading order term, 
as predicted by Theorem [2] and Corollary |3l vanishes if 



T : 



2mT 



/oSi,G(s)ds' 



(71) 



for any integer n. There are three steps in verifying this claim. First, we use Eq. IBHto see that 7g,i(0) — 7g^i(1) 
for the search Hamiltonian. Second, we note that ((5(1)11(1)) = (6(0)11(0)) for the search Hamiltonian. Third and 
last, we substitute these results into Eq. [33]to find that the leading order term in both the upper and lower bounds 
vanish if T is chosen as in Eq. [71] As we only retained the exact expression for the 0{1/T) terms in the error in 
the adiabatic approximation, the higher order terms will not generically be correct. Therefore, we cannot expect our 
error bounds to be asymptotically tight because of the existence of such times. 

On the other hand, if T is not chosen as in Eq. [7T] then we can see that the bounds inferred from Eq. [33] approach 
each other as T increases. This follows because the norms of the derivatives of H{s) are independent of T for the 
search Hamiltonian, and therefore the difference between the two bounds is 0(1/T^). As the upper and lower bounds 
on the error approach each other, we know from the squeeze limit theorem that Eq. [33] and Eq. [33] approach the 
correct expression for the error at all points other than those in Eq. 1711 Thus, our error bounds in Theorem[3]are not 
asymptotically tight in the limit of large T for the Search Hamiltonian; however, they are asymptotically tight for 
any Hamiltonian that is not explicitly a function of T, and for which the leading order contribution does not vanish 
at some T > 0. 
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FIG. 4: This figure shows that our upper and lower bounds given by Theorem [4] for the error in the adiabatic approximation 
converge for a search Hamiltonian with 4 potential marked states, as T — > oo. 

IX. MARZLIN SANDERS COUNTEREXAMPLE 

The Marzlin-Sanders counterexample jl4| is often placed at the center of the ongoing debate regarding the conditions 
that are needed to guarantee the validity of the adiabatic approximation (l6l [l7j . We examine two issues surrounding 
the counterexample, in this section. First, we look at the Marzlin-Sanders counterexample to benchmark our criteria 
for the validity of the adiabatic approximation. We do this to show that our sufficiency conditions prohibit the use 
of the adiabatic approximation for their Hamiltonian. Second, we examine the possibility that Hamiltonians that are 
related to the Marzlin-Sanders counterexample may actually behave adiabatically. In doing so, we provide proofs of 
Corollaries [T] and [H 

The counterexample Hamiltonian is of the form [ij] , 







^ cos(27rs) \ 


Hms{s) = 


UJQ 


sin(27rs) 






^ ) 



2Tr sm{uJosT) 
f 



— sin(27rs) cos{u!osT) 
cos(27rs) cos(a;osT) 
sin(wosT) 



where ax, (Jy and Uz are the Pauli operators. The eigenvalue gap for their Hamiltonian is 

7(s) = 



Wo 



(27rsin(wosT)/r) 

The three most important quantities for determining the error using Eq. [ 

\\dsHMs{s)\\ e 0(1), 

7min e 0(1). 



(72) 



(73) 

(74) 
(75) 
(76) 



By substituting these scalings into our upper bound in Theorem^l we expect the error in the adiabatic approximation 
to be 6(1) and not o(l). This means that our results suggest, but do not imply, that the Marzlin-Sanders counterex- 
ample does not obey the adiabatic approximation. On the other hand, if the second derivative of the Marzlin-Sanders 
Hamiltonian were o{T) instead of ^(T), then the adiabatic approximation would hold for a sufficiently large value of 
T. We noted this fact in Corollary [2l which we prove below. 

Proof of Corollary\^ Since e 0(1) and \\H\\ £ o{T) it follows that F e o{T) and F/T e o(l). We have that 



\H\\T 



r/(7m.nT) 



eo(i) 



(77) 



Therefore Lemma [2] implies that 



11(1 - Fg(1))C/(1,0)|G)|| < Oi + O2 + o(l). 



(78) 
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By applying similar reasoning, we find from Eq. [58] that IIC2II € o(l). As we made no assumptions about ||ij|| in 
the corollary, we cannot use Eq. [Sphere. Instead, we see from Eq . [?T] and Lemma [3] that ||Ci|| € o(l), given that 
\\H\\ e o{T). Therefore if \\His)\\ e o(T) then 

11(1 - PG(l))f/(l,0)|G)|| e 0(1). (79) 

□ 

The above argument shows that some quantum systems can obey the adiabatic approximation if the second derivatives 
of the Hamiltonian diverge sufficiently slowly. We also might imagine that if the derivatives were to diverge even 
slower then it could be possible that the standard criterion for the error in the adiabatic approximation holds for that 
evolution. This claim was made in Corollary [TJ which we prove below using similar reasoning to that used in the 
proof of Corollary [21 

Proof of CorollaryUi Given that \\H{s)\\ e o{VT), Lemma H implies that T/T e o{1/VT). Therefore we have that 



It follows from Eq.[S5]that IIC2II G o(l/T). We then substitute \\H{s)\\ e o{VT) and \\H(s)\\ e o(T) into Eq.EHto 
find that 

\^(mm^7^ ,,(0))T (mm^ 7^ ,,(l))ry 

Hence it follows from Lemma [2] that 

\\H{s)\\ 



(l-PG(l))C/(l,0)|G)||eO max 

\ s=0,l 



(82) 



(min^.7G,i/(s))2r 

completing the proof. □ 



X. CONCLUSIONS 



In this paper, we have shown that the traditional estimate of the error in the adiabatic approximation is valid if 
H{s) is three-times differentiable and the norms of these derivatives are bounded above by a constant for all T > 0. 
We have also provided upper and lower bounds for the error in the adiabatic approximation that are asymptotically 
tight for many Hamiltonians. These bounds are highly relevant to applications in quantum information processing 
because they can be used to provide bounds for the resources needed to perform a quantum adiabatic computation 
using a particular adiabatic path. 

We have also shown that some Hamiltonians that have second derivatives that diverge as T — ^ 00 behave adia- 
batically. Interestingly, the Marzlin-Sanders counterexample Hamiltonian 14] marginally fails the criterion in that 
if the second derivative of their Hamiltonian diverged slightly slower as T — ;> cx), then their proposed Hamiltonian 
would behave adiabatically in that limit. This implies that it is not necessary to forbid Hamiltonians that have second 
derivatives that diverge with T from the class of Hamiltonians that generate adiabatic evolution. 

Our research also leaves open a number of avenues for further work. First, we have noted that for the search 
Hamiltonian, the leading order term in the expansion for the error in the adiabatic approximation vanishes at certain 
times, causing the error to transition from 0(1/T) to 0(l/r^) at those times. This effect could be used to quadratically 
suppress the probability of a particular transition for an adiabatic evolution, and could lead to promising techniques 
for improving adiabatic quantum gates, adiabatic algorithms or quantum control protocols. This effect and its 
implications will be discussed elsewhere. 

Second, a generalization of our bounds for unbounded Hamiltonians that act on an infinite-dimensional Hilbert 
space, or Hamiltonians for which |G(0)) is degenerate, would be an interesting contribution to the understanding of 
adiabatic phenomena. Similarly, our error bounds can be improved in situations where there are near degeneracies 
in the eigenvalues of the Hamiltonians, causing 7niin to be small. In such circumstances it may be more appropriate 
to use integration by parts to approximate all jumps that occur across eigenstates that are separated by a large 
eigenvalue gap, and use a different approximation to handle transitions between nearly degenerate eigenstates. 

Third, our analysis of the search Hamiltonian has shown that the error in the adiabatic approximation appears to 
fall below 2/3 well before the regime where the adiabatic approximation holds. This suggests that bounds for the error 
in the adiabatic approximation are needed in the short-time regime in order to find the time complexity of adiabatic 
algorithms [isj . Such analysis has not yet been performed for adiabatic algorithms, implying that this remains an 
interesting avenue for further research. 
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Appendix A: Bounds for Derivatives of /3^,i, 



We prove our bounds on the derivatives of /? that are presented in Lemma [3l in this appendix. These bounds allow 
us to estimate the contributions to the error in the adiabatic approximation that are incurred due the variation of 
the tunneling amplitude per unit time, (3. Our proof of these bounds is given below. 

Proof of Lemma For convenience we omit the dependence on s of any function below in order to simplify our 
resulting expressions (IE: /S^^^, — /3^,i/(s), 7^,1/ = 7/j,i/(s), • • ■)• Then using this notation, we find by substituting the 



definition of into the LHS of Eg. [39lthat 



can be expressed as 



{i>q\H\iyq^i) {vg\H\iyg_i) 



{Vq\H\Vq^l) %^^^^_^{Vq\H\Vq^l) 



^Vq,Uq-l 



(Al) 



The derivatives of every eigenstate are removed from Eq. lAll by employing a trick involving the resolution of unity. 
For a given eigenstate \v) — \v{s)) we have that 



7V-1 



JV-1 



(A2) 



If {m\iy) = 0, then it follows by differentiating that relationship that (m|i^) = — (to|;>). Similarly, if 1/ ~ m then 
(m|i^) = = — (m|z>), by our choice of eigenbasis. Therefore it follows that 



JV-l 



N-l 



(A3) 



Tn=0 



m=0 



Then by using this property, Eq. IA3[ Eq. IA2I and Eq. [131 the right hand side of Eq. lAll becomes 

(i^q\H\m) (m|ij|i.„_i) , {l^qlHliyq^l) 



E x(^)i^,(5')> E 



lVq,Uq-l 



E 



{Vq\H\m) {m\H\Vq_i) '^^^^^^_^{Vq\H\Vq^l) 



m^Uq-l 



lVq,Vq-l Im.Vq 



lvq,Vq-l 



(A4) 



Eq. Eg. IA4I has now been re-written in a form where no derivatives of \i'q), \i^q-i) nor |m) are present. Furthermore 
the states \vq) and \m) appear as parts of projection operators, which have unit norm. Then by using this property, 
the triangle inequality and lower bounding the eigenvalue gaps in the denominators by 7min we find that Eq. IA4I is 
bounded above by 



2||ij||2 + ||ij||max,,,.,_j7,„,,_J ^ 



7^- 

' mill 



Vq-l,...,lll 



(A5) 



Eq. Eq. [39]thcn follows by using max,,^_,,^_^ I'Tf,,!/,-! I < 2|lij|l. 

By taking the derivative of the inside of the norm in Eq. IA4( and by inserting the resolution of unity to eliminate 
all derivatives of eigenvectors, we find that 



E x(^)k,(s'))/3.,..,-.(s) 
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can be expressed as 



E 

m^i/q 



IVq.Uq-X 



E 

,Vq-ll-^Vq,Vq-l 



lVq,Vq^l 



IVq-.Vq-l 



lVq,Vq-l 



^2 



lVq,Vq-l 



Ivq.Vq^X 



(A6) 



By applying the triangle inequality, using the definition of (3 in Eg. 1131 the fact that for any x and y: ^x,y{s) > 7min, 
< ll^ll/7min> < 2||ij|| and the fact that || i;„Jm(s')) II = 1 for any s and s' , we find 



that Eq. IA6I is bounded above by 



E x(-ik(.'))(|(e 



/^i^q , m /^m , i/g _ 1 O^i/g , m V ^ I^Uq^ml^m^Vq — x^ra.Uq-x 



lfq,Vq-l 



E 

m^Vq 



lVq,yq-l 



E/^M / t\\ ( "f'^q,'^q-l(^'^q-Mq-l ~^ ivq.Vq-lPuq,Vq-l 



Vq,Vq-l,...,Vl 



\H\\ , 4||i7|| 



lVq,Vq-l 



7'- 

' mm 



^min Tniin 



(A7) 



Our final result in Eq. HOlis found by applying the triangle inequality to this, using Eg. 1391 and the following result 
from [H, 



< 2||H|| 



8||g|| 



(A8) 



We then find Eq. 1301 by using the triangle inequality and equations Eq. IA8I and Eq. to upper bound Eq. IA61 □ 



Appendix B: Proof of Lemma |4] 



We present a proof of the upper bounds that are presented in Lemma 2] for the contribution to the error that results 
from the sum over all one- and two-jump paths. Our proof technique is similar to that used in the proof of Lemma [2] 
The most significant difference is that we have to apply integration by parts twice to approximate the sums of all one- 
and two-jump terms, which we denote Ci and C2. 

Proof 0/ Lemma 

The contribution of the sum of all one-jump paths to the error in the adiabatic approximation is 



JV-l 1 

Ci = E 1/^(1)) / /3M,G(s)e- 
u=l -^0 



(Bl) 



We find the following upper bound on || Ci || by applying integration by parts twice and the triangle inequality on Eq. lBll 



l^^ill < 



Ei'^iW) 



El'^i(l)) 



d 



ds -ilG,^i{s)T \ds -ilG,uAs)T 



(B2) 
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Using our bounds for the derivatives of /3 in Eq. |39]and Eq. |40l the bound |7iyi,iy| < 2||ij|l and the upper bound for 
I7I in Eq. [M] wc find that 



l^^ill < 



Rn 



where 



2||g|| + ||g|| , 20||g||||g||+12||gf , 88\\Hf 
■K-o Q — — 1 A — — H 



73. T2 

' mm 



T2 

/ mm 



7^. r2 

/ mm 



(B3) 



(B4) 



Now we win apply similar reasoning to approximate C2, which is the contribution of all non-adiabatic two-jump 
paths to the error in the adiabatic approximation. This error can be expressed as, 



We use integration by parts once on the outer-most integral and the fundamental theorem of calculus to find 



(B5) 



^ ( d /3i.2,i'i(si) 



(B6) 



We must have that E^^{s) 7^ Eq{s) to use integration by parts again on this expression, which is ensured by our 
assumption that the initial state is non-degenerate. Using integration by parts again to estimate the inner-most 
integrals in Eq. IB6I we have that 



E 1^2(1)) 



I ds2 ~i"fG,ui{s2)T j 



S2 



\ dsi -i-iut.v2{si)T 



\ -ilG.uAs2)T 



, ( A^-£^)e-/o-7o,.(OdCTd, L., 

\ds2~l"fG,yAS2)T 



/3.„.,(si)/3,„G(si)e-'/o^7G..,(«)+7.„.,dCT 



(si)(7G 

^ ( _d /3^2..iyiisi)P:^l,Gi^l) 

\ OSi —71/1,1^2 (si)(7G ,1/1 (si) + '7^1,1^2 



(B7) 



Then using the triangle inequality, our upper bounds for the derivatives of (3 in Eq. [39] and Eq. 1401 the bound 
\ii^2,i^\ ^ and the upper bound for I7I in Eg. IA8I we find by substitution that 



1^ I, ^ ||g|P+4||H||2 + 5||g||||g|| ^ 12||g||||H||2 + 30||gf ^ 
\'^2\\ S -A — — \ s — — H 



74. T2 

'mm 



rp2 

I mm 



76. r2 

' mm 



(B8) 
□ 
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